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Magnetic fields permeate the entire solar atmosphere weaving an extremely com¬ 
plex pattern on both local and global scales. In order to understand the nature of 
this tangled web of magnetic fields, its magnetic skeleton, which forms the bound¬ 
aries between topologically-distinct flux domains, may be determined. The magnetic 
skeleton consists of null points, separatrix surfaces, spines and separators. The skele¬ 
ton is often used to clearly visualise key elements of the magnetic configuration, 
but parts of the skeleton are also locations where currents and waves may collect 
and dissipate. 

In this review, the nature of the magnetic skeleton on both global and local 
scales, over solar cycle time scales, is explained. The behaviour of wave pulses in 
the vicinity of both nulls and separators are discussed and so too is the formation of 
current layers and reconnection at the same features. Each of these processes leads 
to heating of the solar atmosphere, but collectively do they provide enough heat, 
spread over a wide enough area, to explain the energy losses throughout the solar 
atmosphere? Here, we consider this question for the three different solar regions: 
active regions, open-field regions and the quiet-Sun. 

We find that the heating of active regions and open-field regions are highly 
unlikely to be due to reconnection or wave dissipation at topological features, but 
it is possible that these may play a role in the heating of the quiet-Sun. In active 
regions, the absence of a complex topology may play an important role in allowing 
large energies to build up and then, subsequently, be explosively released in the 
form of a solar flare. Additionally, knowledge of the intricate boundaries of open- 
field regions (which the magnetic skeleton provides) could be very important in 
determining the main acceleration mechanism(s) of the solar wind. 
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1. Introduction 


The question of how the solar corona (or indeed any stellar object with a hot 


corona) may be heated has been considered for many decades (see, e.g. Parnell & 


De Moortel, 2012 for a recent review), but still remains unanswered. One significant 


advance has been the recognition that explaining simply how the corona alone is 
heated is insufficient; the real question is how does the whole solar atmospheric 
system (the photosphere, chromosphere, transition region and corona) interact and 
interlink in order to sustain a hot corona? Throughout this volume there are many 
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articles that look to address different aspects of this problem. In this article, we 
focus on magnetic topology and the role it plays (if any) in the heating of the solar 
atmosphere. 

Magnetic topology incorporates all properties of a magnetic field that are pre¬ 
served by ideal displacements. For instance, examples of topological features are: 
the linkage and knottedness of field lines, null points, their associated separatrix 
surfaces and spines, as well as separators (see Sec.|^. The topology of a magnetic 
field is not changed by stretching it; it can only change through the process of 
reconnection. Indeed, a change in topology implies reconnection. 

Almost a decade ago a review article was published on aspects of magnetic 
topology (Longcope 2005). Since then, there have been many developments on 


two fronts: (i) our understanding of the nature of the topology of both global and 
local coronal magnetic fields and (ii) the role particular topological features play in 
processes such as magnetic reconnection or wave dissipation. Before reviewing these 
developments, the basic elements of 3D magnetic topology are outlined in Sec. 
Then, in Sec. we investigate the nature of the magnetic skeleton of the global 
coronal magnetic field, and how it varies over the solar cycle, as well as consider how 
the magnetic skeleton may evolve under magnetohydrodynamic (MHD) conditions. 

The energy required to heat the solar atmosphere is injected through the pho¬ 
tosphere from the convection zone below. This energy arrives in two main forms: 
rapid motions which incite waves to travel up into the atmosphere or slow motions 
that cause a gradual stressing of the magnetic field. In Sec.|^ we discuss the impact 
these different driving motions have on the elements of the magnetic skeleton. Fi¬ 
nally, in Sec.[^ we discuss what the consequences of these results are for heating the 
solar atmosphere and, thus, address the question: is magnetic topology important 
for coronal heating? 


2. Basic Elements of the 3D Magnetic Skeleton 

In order to determine the complexity of a magnetic field in 3D, it is beneficial to 
determine its magnetic skeleton; this comprises of several features. Magnetic null 
points are locations where all three components of the magnetic field equal zero. 
An infinite number of field lines extend from(/to) a null point forming a surface, 
known as a separatrix surface^ and a pair of lines extend into(/out from) the null, 
known as spines. If the field lines in the separatrix surface are directed out of the 
null, then the spine field lines will be directed inward, and the null is known as 
a ‘positive null’. If these field line orientations are reversed, the null is said to be 
‘negative’. For true 3D magnetic null points no other scenarios are possible, as the 


magnetic field must be divergence-free (e.g. Parnell et al. 

1996 

. Both the local 

magnetic field near the null and field lines in the separatrix surfaces can take on a 

range of geometries; these were categorised by 

Parnell et al. 

(1996 

). 


Special field lines, called separators, may link connecting pairs of null points. 
These can be formed in several different ways, but the only generic type of separator 
is that formed by the intersection of two separatrix surfaces from opposite-polarity 
null points. Fig. illustrates each of these features in a simple magnetic field in¬ 
volving two null points (one positive, the other negative), their associated separatrix 
surfaces and spines, and a separator linking the two nulls. 
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Figure 1: A comparison of the same magnetic field, but viewed in two difi'erent 
ways by: (a) drawing the topological features in the field and (b) drawing sample 
field lines in four different colours. In (a) the topological features include: a pair 
of oppositely-signed null points (pink/blue spheres), their associated spines (thick 
lines) and separatrix surfaces (pink/blue field lines and shaded surfaces). These two 
separatrix surfaces intersect along a line called a separator (green line) that links 
the two oppositely-signed null points. 



Figure 2: Example of a coronal loop structure which has been interpreted as the 
magnetic field about a (2D) X-type null configuration, but it could equally be 
the magnetic field about a 3D separator (c.f. Fig. [^). This 17lA image, taken by 
TRACE, shows coronal loops from the active regions NOAA 9149 (above) and 9147 
(below), taken at 10:17UT on 4**' September 2000. 


Fig.[T)D illustrates the same magnetic configuration as that shown in Fig. 
but only for several field lines which do not approach the separator. It is interesting 
to note that these field lines map out a feature, which, if viewed from above, would 
look like the field about a 2D X-point. There are a number of observations of coronal 
loops that may be/are interpreted as X-point like; an early example from TRACE 
is shown in Eig. In general, similar coronal loop formations have been attributed 
to there being a null point at the ‘X’ (e.g. Narukage et al. 2014 Freed et al.. 2015), 


however, it is possible that a separator may be present instead. 

The local separator magnetic field in planes perpendicular to separators is not 
always X-type, but may be 0-type, if there is a sufficiently large component of 
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current parallel to the separator (Parnell et al ., 2010a). In such cases, any magnetic 


field lines that have undergone separator reconnection (thus having presumably 
been heated) and so are potentially visible in coronal images, would not show an 
0-type form. This is because these field lines are actually helical in nature (with 
a much stronger magnetic field component along the separator than the 0-type 
component perpendicular to it) and, therefore, twist about the separator. So in these 
cases the expected signature would be a narrow, long, twisted loop-like structure. 

Null points in the corona are, of course, 3D not 2D, but in many instances they 
may appear 2D, if the field lines in the separatrix surface are strongly aligned along 
one line (for example, if the local field about the null has a major fan eigenvalue, as 


explained in Parnell et al. 19961. More importantly, null points are local magnetic 


structures that are identified by determining where B = 0. So, by definition, they 
are single points. Separators, on the other hand, are global magnetic structures and 
cannot be identified directly by local measurements of the magnetic field. They form 
lines and, thus, have a far greater spatial extent than null points. In particular, as 
shown in Sec. there are many large-scale separators that arch through the solar 
atmosphere achieving lengths that are longer than a solar radius. 


In the recent work of Freed et al. (20151, the location of coronal null points 


identified from low-resolution potential field source surface (PFSS) extrapolations, 
are compared with AIA observations of the Sun’s coronal field. In about 31% of the 


cases, Freed et al. (2015) found that the AIA observations showed structures, such 
as X-type loop patterns, that could be interpreted as the configuration associated 
with a null point. However, since separators are field lines that generally extend 
from or to a null point, it is not unreasonable to imagine that in a proportion of the 
above cases, the coronal signature may reveal the location of a separator, rather 
than a null point. It would be interesting to redo this analysis with the locations of 
the separators also identified. 


3. 3D Coronal Topologies: Global &: Local 

All topological features are fundamentally associated with magnetic null points[t] 
The numbers of coronal null points that exist in potential magnetic field configura¬ 
tions, created by direct extrapolation from observed quiet-Sun regions, have been 


determined exactly by locating all the nulls (e.g., Regnier et al. 2008) or estimated 


using a spectral method (e.g., Longcope & Parnell 2009 Longcope et al. 2009). 


Recently, the nulls that occur in PFSS models of the global coronal magnetic field 


(Cook et al. 2009 

Platten et al. 

2014 

Edwards 

2014 

Freed et al. 

2015 

Edwards 

& Parnell 

2015) have been counted. 


Studies, whose focus is the determination of the number of separators or even 


the complete magnetic skeleton, are uncommon. Close et al. (2004) investigated the 


t In dom ains where the magnetic field is not closed, separators may instead connect to bald 
patches (e.g. [Titov et 11 /"] |1993| |2011[ |. This type of behaviour occurs where there is a null point 
lying outside the domain. 
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Figure 3: The global coronal topology of two PFSS magnetic fields constructed 
from high-resolution SOLIS synoptic magnetograms taken at (a) solar minimum 
(CR2083) and at (b) solar maximum (CR2130). The PFSS models use spherical 
harmonics (with a maximum harmonic number of I = 301) to extrapolate from the 
magnetogram data. Topological features marked are: null points (red/blue spheres), 
spines (orange/purple lines) and separatrix surfaces (pink/blue thin lines). Where 
the separatrix surfaces intersect the source surface, at i? = 2.5Rq, there are thick 
pink/blue lines. The thick green lines on the source surface indicate the base of 
the HCS, whilst the thin green lines extending down from these map out the HCS 
curtains, which divide open/closed field from the Sun. Cuts, at constant radii of (c) 
& (d) i? = 2.5i?0 and (e) & (f) i? = l.OOSi?© through the topological features, are 
plotted for the solar minimum case on the left and the solar maximum case on the 
right. Overplotted on (e) and (f) are all the null points found in each case. 
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separators found in regions of quiet-Sun, but the most comprehensive survey to 
date of magnetic skeletons is by Flatten et al. (2014). 


(a) Topological Structures of Global Coronal Fields 


Flatten et al. (2014) determined the numbers of coronal null points and sepa¬ 


rators in 496 global coronal potential fields, taken once every Carrington rotation 
(CR), spanning three solar cycles. The potential fields were constructed from spher¬ 
ical harmonics (with a maximum harmonic number I = 81) using a FFSS model, 
where synoptic magnetograms, taken at Kitt Feak by the NSO Vacuum Telescope 
and SOLIS (low resolution), were used to define the magnetic field on the solar 
surface and the source surface was taken at R = 2.5 Rq. Two example magnetic 
skeletons of global fields constructed from high-resolution SOLIS magnetographs 
are shown in Fig. 3] _ 

The study by Flatten et al. (2014) revealed a number of interesting features. 


(i) Null points above active regions occur at much higher heights than quiet-Sun 
nulls, (ii) The number of coronal nulls varies out-of-phase with the solar cycle. 
In the two examples shown in Fig. there are 1964 nulls in the solar-minimum 
case, CR2083, and almost half that number, 1131 nulls, in the solar-maximum 


case, CR2130. This result is the opposite of that found by Cook et al. (2009) and 


Freed et al. 

(2015 

) who find that null-point numbers vary in phase with the solar 

cycle. However, neither of these studies consider small-scale magnetic features. In 

Cook et al. 

O 

o 

, simulated (rather than observed) magnetograms are used from 


a flux transport model in which only active regions (i.e. large-scale structures) are 
emerged into an otherwise smooth magnetic field. In Freed et al. ( 2015| ), the global 
extrapolations are from observed magnetograms and have a maximum harmonic 
number of just I = 30, so small-scale structures are missed. This means these two 
studies can only reliably find nulls at high altitudes, which preferentially occur 
above active regions, (iii) More null points are found during solar minima with 
weak polar fields than minima with strong polar fields. 

During solar minima, there are large expanses of quiet-Sun photospheric field 
that are filled with small-scale magnetic features of both polarities. These produce 
a complex tangled network of closed coronal field which contains many null points 
(see Fig. These quiet-Sun fields are surrounded by the dipolar field originating 
from the Sun’s polar regions. If the polar fields are strong then the dipole field will 
dominate, providing a strong confinement of the quiet-Sun fields and thus limiting 
quiet-Sun nulls to low heights. On the other hand, if the polar fields are weak, then 
the dipolar field does not constrain the quiet-Sun field allowing the nulls within it 
to reside at higher heights above the solar surface. This indicates that during weak 
polar-field minima (such as the minimum period between cycles 23 and 24) nulls 
are likely to occur at higher altitudes. Nulls at higher altitudes are more reliably 


identified, hence, in the work of Flatten et al. (2014) more nulls are found during 


the weakest polar-field minimum investigated. 


The separators found in global FFSS models come in two forms (Flatten et al. 


2014): those that connect pairs of coronal nulls (called null-null separators) and 


those that connect a coronal null to the ring of nulls that forms the base of the 
heliospheric current sheet (HCS) on the outer boundary (source surface) of the 
model. Figs. H and ii show cuts at the source surface, through the magnetic 
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skeleton of the two global fields shown in Figs.|^ and[^. The HCS null line is the 
thick green line in these cuts. There are more null-null separators found during solar 
minimum than during solar maximum. For instance, 1997 separators are found in 
CR2083 (solar minimum) of which 1946 are null-null and 51 are null-HCS and just 
808 separators, in CR2130 (solar maximum) of which 765 are null-null and 43 null- 
HCS. The number of null-HCS separators is not strongly dependent on the solar 
cycle. 

Finally, we note that not only do separatrix surfaces emanate from individual 
coronal nulls, but they also extend out from the HCS null line (green lines in 
andj3|) separating regions of open field from regions of closed field 


Figs. J^,[ 

e.g. 


(e.g. [Titov et al. 2011). At solar minimum, the HCS null line is approximately 


equatorial, however, at solar maximum, it is highly warped and may cross the poles 
(Fig. It!) ■ Furthermore, it can split into two or more disconnected loops (e.g. |Wang| 


et al. 2014[ Edwards et al. 2015). 


The HCS null line arises because the magnetic field outside the source surface 
is assumed to be purely radial: where the radial field changes direction a null line 
forms. In reality, the magnetic field of the solar atmosphere does not all become 
purely radial at a single radius above the surface and so a null line will never occur in 
practice. Instead, additional coronal nulls will occur with their associated separatrix 
surfaces. This means that in reality there will be no null-HCS separators, just 
additional null-null separators connecting nulls low down in the solar atmosphere to 
nulls lying above 2.5i?Q. It is not immediately clear how many null-HCS separators 
currently found will become null-null separators, as opposed to simply disappearing, 
nor is it clear whether additional null-null separators will be created between nulls 
above 2.5Rq and nulls below 2.5Rq that currently do not have null-HCS separators. 
Further work is required to understand the nature of far coronal magnetic field. 


(b) Topological Structures in MHD Experiments 


Studies of magnetic topology and null point numbers typically focus on potential 
magnetic fields. However, there have been a small number of studies investigating 
the magnetic skeletons of magnetic fields created in highly-dynamic numerical MHD 


experiments ( 

Haynes et al. 

2007a 

Maclean et al. 

2009 

Parnell et al. 

20106 

Mac- 

Taggart & Haynes 

2014). By highly dynamic, we mean that the magnetic fields 


undergo significant changes, rather than they evolve at a particularly fast rate. 


The experiment considered by Maclean et al. (2009) and Parnell et al. (20106) 


involved the emergence of a twisted flux rope from below the solar surface up into 
an overlying horizontal coronal magnetic field angled at 135° to the axis of the flux 


tube (Galsgaard et al. 2007). In the experiment, the centre of the flux tube was 


made buoyant such that it rose up and interacted with the overlying field. Initially 
there were no nulls within the system, but the interaction between the flux tube 
and coronal field led to the creation of two null clusters, each containing multiple 


null points, either side of the emerged flux tube (Maclean et al. 2009). The null 


points within each null cluster were linked together by single short separators, like 
beads on a string. All but four nulls were short lived, with the associated cluster 
separators being equally short lived as the nulls they connect. The long-lived nulls 
(two of the same sign inside each null cluster) were connected by a multitude of 
long separators (up to several hundred) that link from one null cluster to the other 
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arching up over the emerged flux tube (Fig.|^). Each separator lies at the interface 
between four topological domains containing: overlying field, flux-tube held and 
held that connects from the hux tube to the overlying held and vice-versa. These 
long separators, called intercluster separators, typically survived much longer than 
the short-lived nulls and cluster separators, but only a few had lifetimes comparable 
to that of the long lived nulls. 



(C) “ ” ” (d) 

Figure 4: (a) Field lines and elements of the magnetic skeleton in a frame from a 
hux-emergence experiment ( [Parnell et al. 201061: pink/blue lines forming the hux 
rope with some also connecting to the overlying held. The null points (red/blue 
spheres) and separators (coloured according to the parallel electric held (low - 
black, high - red) along them) are also plotted, (b) Connectivity of the magnetic 
held threading a vertical plane perpendicular to the hux rope half-way along its 
length: overlying held (red), hux rope (blue), hux rope-overlying (yellow), overlying 
to hux rope (cyan), (c) Integral of the electric held parallel to the magnetic held 
threading the same plane (white - no reconnection and red/purple - greatest amount 
of reconnection). (d) Regions of high Q mapped onto the same plane (white - Q = 
0 and blue - high Q). Figure based on images from Restante (20111. 


Restante ([2011 ) undertook a thorough analysis of the topology, quasi-separatrix 


layers (QSLs, see e.g. 

Priest & Demoulin 

1995 

Demoulin et al. 

1996 

Titov et al. 

20021 Titovl 2007 

1, and the primary sites of reconnection in 

this 3D numerical 


experiment. QSLs are geometrical rather than topological features of the magnetic 
held. Neighbouring held lines within a QSL will have foot points that are far apart. 
The extent of the separation of the foot points is measured by the squashing factor, 
Q, which is designed, such that in a given plane, each end of the same held line has 


the same value of Q (Titov et al. 20021. Regions of high Q have been associated 


Titov et al. 

21 mi;; 

Galsgaard et al. 

2003 

Aulanier et al. 


2005, 20061. 
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Several field lines (seen in pink/blue in Fig. |^) form a twisted flux rope that, 
at the time of this frame, had emerged and partially reconnected with the overlying 
field. This is illustrated in Fig. by the few flux-tube field lines that are open and 
extend up into the solar atmosphere above. Also seen in Fig.|^ are the null points 
(red/blue spheres) in the two null clusters and the separators connecting the nulls 
within the same null cluster (thick black lines) and connecting nulls in different null 
clusters (thick multicoloured lines). 

To compare the magnetic skeleton, QSLs and sites of reconnection specific quan¬ 
tities were calculated on the plane that cuts through the flux tube at right angles, 
half-way along its length: this plane was chosen because it intersects all the field 
lines in the volume. Separatrix surfaces denote the boundaries between regions of 
different connectivity. In Fig.[^, the intersection points of field lines with this plane 
are coloured one of four colours according to their connectivity: overlying field (red), 
flux rope (blue), flux rope-overlying (yellow), over lying-flux rope (cyan). The sites 
of reconnection (Fig. [^) were identified by plotting contours of the integral along 
a field line of the electric field parallel to the field line {J E^^dl - a, non-zero value of 
this integral is a necessary and sufficient condition for 3D reconnection 


Schindler 


et al.. 1988). Finally, in Fig. |4p, the value of Q on each field line is indicated, to 
identify sites where field lines, that start out close to one another, diverge away 
significantly from each other. Note, to determine Q, both ends of the field line must 
intersect the photosphere of the model. So, Q is undefined for all overlying field 
lines, or field lines connected, at one end or the other, to the overlying field. 

Comparing Figs. @3# and[^, it is clear that the sites of reconnection coincide 
well with boundaries between different connectivities, i.e., they coincide with topo¬ 
logical features. However, even though there are regions of high Q in these locations, 
there are also several other large regions of high Q. In these additional high-Q re¬ 
gions, there is no reconnection. By a careful analysis of the magnetic field in these 
regions, Restante (2011) found that they coincided with significant changes in field 
line geometry. Therefore, it is important to be cautious about the interpretation of 
QSLs. 

Not all QSLs are sites of reconnection. But it is also the case that elements in the 
magnetic skeleton are not always sites of reconnection. In order for reconnection 
to occur, both a favourable magnetic field configuration and favourable plasma 
conditions are required to create the non-zero parallel electric field that is essential 
for 3D reconnection. A key question, therefore, is: Can magnetic field conditions 
favourable for reconnection arise in the absence of either QSLs or elements of the 
magnetic skeleton? 

Firstly, from a practical point of view, in order to determine Q on a field line, 
both ends must pass through the surface upon which the separation of originally 
nearby footprints is measured. This is not always the case, as seen in Restante] 
(2011) where Q was only defined for field lines that crossed the photosphere at both 
ends, i.e., Q was found for flux rope field lines, but was undefined for flux rope- 
overlying, overlying-flux rope or purely overlying field lines. Some of these field lines, 
for which Q could not be defined, are associated with reconnection (as evidenced 
by the wings in Fig. |^). In this experiment, this reconnection is associated with 
topological features, but it is possible that, in a different situation where there are 
no topological features, similar field lines where Q can not be defined may undergo 
reconnection. 
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Secondly, there are many reconnection experiments, which have not been anal¬ 
ysed to determine one or all of the QSLs, the magnetic skeleton or the / E\\dl. 
Therefore, it is not possible to say where exactly the reconnection is occurring and 
what it is associated with. In one reconnection experiment by |Wilmot-Smith et al.\ 
(2009), a very detailed analysis was undertaken. They found that although Q and 
f Ej I dl both had a similar form, the regions of highest Q and highest J Ej^dl did not 
coincide. Further studies of reconnection experiments undertaking similar detailed 
analysis are required before the exact relationship between Q and f E^dl may be 
properly understood. 


4. Responses to Photospheric Driving 

It is well known that the source of energy that heats the solar atmosphere and 
powers events, large and small, from solar flares and coronal mass ejections down 
to X-ray bright points and nanoflares, comes from convective motions below the 
solar surface, which either inject energy during the process of flux emergence, or 
inject energy by driving photospheric magnetic footpoints. The footpoints may be 
driven at two different rates. Motions that are faster than the Alfven speed generate 
waves which may propagate into the solar atmosphere, whereas motions that are 
slower than the Alfven speed cause a stressing of the magnetic held, throughout the 
atmosphere, that can lead to the generation of current layers. Below, we focus on 
what happens at elements of the magnetic skeleton in response to these motions. 


(a) Behaviour of Wave Pulses at Topological Peatures 


In recent years, the advent of telescopes with both high spatial and temporal res¬ 
olution, as well as high sensitivity, such as SoHO/EIT, Hinode/XRT, Hinode/EIS, 
SDO/AIA and CoMP, have enabled a vast number of wave observations to be made. 
What happens to waves when they come across features such as nulls, separators 
or separatrix surfaces? What happens to these topological features when they are 
hit by waves? 

Firstly, we note that the waves which travel across the magnetic field are fast- 
magnetoacoustic waves. This type of wave travels at a speed, c/ = where 

VA = B /is the Alfven speed {B is the magnitude of the magnetic field, p is 
the plasma density and p, is the permeability of free space) and Cs = \/pJp is the 
sound speed {p is the plasma pressure). 

The solar corona is a low plasma-beta environment, where the magnetic pressure 
dominates over the plasma pressure, so va^ Cg', thus, Cf ^ v a- At a coronal null 
point, 5 = 0, therefore va = 0 suggesting that a fast-mode wave will be unable to 
propagate across a null. 


McLaughlin & Hood (2004) explored the behaviour of fast-wave pulses in the 


vicinity of a 2D null, in a zero-beta plasma, and showed that the wave essentially 
becomes trapped. However, the speed of the wave also depends on the sound speed. 
At a coronal null point Cf = Cg, which, although small, is non zero and, hence, 
a fast-mode wave never actually gets trapped at the null. Instead, the wave un- 


Hood 

(2006 

) ano 

Thurgood & McLaughlin 

(2012 

1 and in the Hall MHD regime by 

Threlfall et al. (1 

10121. In all cases, the trapping or mode coupling of the fast-mode 
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wave gives rise to localised oscillatory reconnection at the null itself, as well as 
shocks. Both of these cause heating along the field lines that pass close to the null. 

Recently, [Thurgood fc McLaughlin ( 2013&|a ), have studied the behaviour of 
Alfven waves (which travel at the Alfven speed parallel to the field), in the vicinity 
of both 2D and 3D nulls. Their work suggests that these waves may become phase- 
mixed near null points, again causing heating. 

As discussed above, both Alfven and fast-magnetoacoustic waves can be affected 
by the presence of null points (for a review see McLaughlin et al ., 2011). But what 
happens to waves as they approach separators or separatrix surfaces? This question 
has never been addressed before, but since the magnetic field does not reduce in 
strength near a separator or separatrix surface, it is quite possible that these waves 
simply cross these features without any change in behaviour. It would still, though, 
be interesting to investigate what, if any, modification in behaviour is found across 
separators and separatrix surfaces to confirm or deny the above hypothesis. 

In this section, we have so far addressed the behaviour of just single wave pulses 
interacting with null points. However, analytical and numerical modelling has shown 
us that single pulses in MHD (Hood et al. 20021 and in Hall MHD (Threlfall, 20121 
are often limiting cases of full wave trains. So these works may be seen, in some 
sense, as special cases of the situations we discuss below that undergo continuous 
driving on the boundary. 


(6) Effects of Boundary Driving on Topologieal Features 


The work discussed above considers what happens when a single wave pulse 
arrives at a null. How do topological features respond to continuous stressing as 
opposed to a single wave? In situations where the local field about the null is driven 
slowly away from its initial potential state, waves are created that lead to the 
collapse of the null, creating a localised current layer and triggering reconnection. 
An enormous body of work exists studying reconnection at 2D null points triggered 
by boundary driving (e.g. see reviews by [Priest fc Forbes 2000 Biskamp 2000 and 
enclosed references). Similarly, reconnection at 3D null points has, in recent years, 
received quite a bit of attention. 

Due to the additional complexity that the third dimension brings, in 3D, the 
specifics of the boundary perturbation can cause a variety of different current accu¬ 


mulations, that are associated with different types of reconnection (Priest & Pontin 


2009). In particular, a 3D potential null point’s spine is orthogonal to its separatrix 


surface. If the angle between the separatrix surface and spine is decreased then this 
can precipitate a collapse of the null, in the same manner as a 2D null collapses, 
with a component of current created perpendicular to the plane of the collapse. 
This causes spine-fan reconnection. If, instead, the separatrix surface is perturbed 
in such a way that the spine and separatrix surface remain perpendicular, but the 
field about the spine or about the plane of the separatrix surface is twisted, then 
instead a current component parallel to the spine is formed. Here, either torsional 
spine or torsional fan reconnection occurs. The reconnection triggered is different in 
all three cases. A comprehensive review of 3D null-point reconnection can be found 


Priest & Pontin (2009) and Pontin et al. (2013). Examples of solar scenarios 


in which reconnection has been triggered at null points by boundary driving in¬ 


cludes work by Aulanier et al. (2000); Pariat et al. (2009); Masson et al. (2009) and 
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Masson et al. (2014), but these cases are all aimed at modelling coronal mass ejec¬ 


tions (CMEs) or eruptive events, such as solar flares, rather than studying coronal 
heating. 

The response of magnetic separators and separatrix surfaces to similar boundary 
driving has received less attention than that for null points (probably because sep¬ 
arators and separatrix surfaces are harder to find). However, Haynes et al. (2007&); 


Parnell et al. (2008 2010a| and Parnell et al. (2010&) study the evolution of the 


whole magnetic skeleton in driven MHD experiments. Haynes et al. (200761 showed 
that boundary driving leads to the build up of currents about the separator and 
that the resulting reconnection triggered in the current layer naturally leads to 
the creation of multiple separators connecting the same pair of nulls. |Parnell et al.\ 
(20106) also found this behaviour in a completely different system. In both cases, 
the increase in complexity was required to enable a rapid rate of reconnection, and 
thus, was also connected with the most intense period of heating. As the recon¬ 
nection rate decreases, so too does the complexity of the magnetic skeleton. This 
enhancement in complexity/mixing appears to be a generic feature of fast 3D recon¬ 
nection and has been found in other MHD reconnection experiments, in the absence 
of null points (Pontin et ail 2011), as well as at null points (Wyper & Pontin 20141 


and also in a kinetic reconnection experiment (Daughton et al. 2014). 


(c) Energy Partitioning due to Reconnection at Topological Features 


In the above studies, reconnection was triggered by boundary driving. Here, we 
consider the studies that have taken a different approach and start instead from 
force-free and non-force-free MHS equilibria involving current layers where recon¬ 
nection is initiated by micro-instabilities (in MHD experiments, an anomalous re¬ 
connection is applied to mimic the triggering of reconnection by micro-instabilities). 

Large Lorentz forces are created within the strong current layers found in non¬ 
force-free equilibria. These are counter-balanced by gradients in pressure, providing 
overall force balance |t] 

One interesting aspect of these reconnection studies is the ability to follow, in 
great detail, the energy partitioning and transport resulting from the reconnection. 


Longcope & Priest (2007) and Longcope & Tarr (2012) considered what happens 


following reconnection in a force-free equilibria with a current sheet located at a 2D 
null point. They analytically determined that the proportion of magnetic energy 
going directly into Ohmic heating was much smaller than that converted into kinetic 
energy and ultimately released by viscous dissipation. Thus, they found that 2D 
magnetic reconnection (in an initially force-free system) could cause significant 
heating over a wide region, potentially far from the null point. 

On the other-hand, investigations of the energy partitioning resulting from re¬ 
connection in non-force-free MHS equilibria with a current sheet at a 2D null point 


by Fuentes-Fernandez fc Parnell|| (2012) and Fuentes-Fernandez & Parnell (2013) 


find that the kinetic energy created in the system is an order of magnitude less 


t Null points and separators undergo an infinite time collapse, therefore, it is not possible to 
reach a perfect force balance about such features. Instead, a state is reached in which the total 
force is zero everywhere, except at the nulls and separators themselves, where it is very small (e.g. 
[Fuentes-Fernandez &; Parnellj [2012[ [^13[ [Stevenson et q/.[[^15[ |. 
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than the changes in magnetic and internal energy. In these studies, Ohmic heating, 
directly in the vicinity of the null point, dominates over all other forms of heating. 

Both high and low plasma-beta cases have been considerecjf] with more energy 
going into the waves in the low-beta case, than the high-beta case. These waves are 
created by the sudden loss of force balance and emanate out from the edge of the 
diffusion region. They are magnetoacoustic in nature and travel at the fast-mode 
speed everywhere except, along the separatrices (which are field lines), where they 
travel at the slow-mode speed. 



Figure 5: (a) Contour plot of the current in a non-force-free, MHS equilibria about 
a 2D null point, (b) Perturbation in pressure indicating the waves and flows created 


by the localised reconnection at this 2D null (created from Fuentes-Fernandez & 


Parnell 2013). (c) The magnetic skeleton of a non-force-free 3D MHS equilibria 


involving a separator current layer (indicated by an isosurface), (d) Perturbations in 
pressure in the xy-plane half way along the separator indicating the waves and flows 


created by reconnection in the current layer (created from Stevenson & Parnell 


2015). 


Stevenson & Parnell (2015) considered a system involving a non-force-free 3D 


MHS equilibria in which a separator current layer exists with current parallel to 
the separator. Only high-beta cases have been considered thus far; in all cases the 
waves and flows generated by the reconnection carry just a fraction of the magnetic 
energy converted during the reconnection away from the reconnection site. It would 
be interesting to investigate a low-beta case of separator reconnection, where, there 
would be no issues about the plasma beta becoming infinite, instead remaining low 
throughout the reconnection site and surrounding region. 

t Exactly at the magnetic null point the plasma beta will be infinite no matter what value the 
plasma pressure is there. 
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The reconnection proceeds in essentially the same way as it does at a 2D null 
point. Rapid reconnection occurs within the separator current layer, converting 
most of the released magnetic energy directly into internal energy via Ohmic heating 
at the separator. This causes a loss of force balance which generates waves that are 
launched from the edges of the diffusion region. In planes perpendicular to the 
separator, these waves behave in a very similar fashion to that seen in the 2D null 
point case (Fig. [^. 

In the wake of these waves, the system tries to regain force balance and so 
flows are generated. These flows drive slow, steady reconnection at the separator, 
which lasts much longer than the initial fast-reconnection stage. From [Stevenson*^ 


Parnell (2015), it is clear that the separator reconnection does not occur at the 3D 


null points that lie at either end of the separator, but occurs at some position along 
its length. The initial seed location of the reconnection is likely to be where the 
separator current reaches a peak. This could be anywhere along the length of the 
separator, depending on the nature of the surrounding magnetic field and plasma. 

The partitioning of magnetic energy released by reconnection is still an unan¬ 
swered question, even in the relatively simple MHD regime where the magnetic 
energy can only become kinetic or internal energy. In reality, it could also be con¬ 
verted into non-thermal energy producing accelerated particles and thus even less 
will then go into internal or kinetic energy. In order to explain the solar atmospheric 
heating problem, it is essential that the energy partitioning due to reconnection is 
properly understood, whether the reconnection occurs at topological features or in 
their absence. 


5. Discussion &: Conclusions 


In this paper, we have given a brief review of both the nature of the topology 
throughout the corona, as well as discussing the different energy release mecha¬ 
nisms associated with the elements of the magnetic skeleton. Here, we consider the 
question posed in the title: is magnetic topology important for coronal heating? 
The answer to this is not straightforward and depends on a number of things, (i) 
how the question is interpreted, (ii) for what location in the solar atmosphere, and 
for what period during the solar cycle, the question is posed. 

Let us first consider whether the trapping, coupling and dissipation of waves at 
null points or reconnection at null points, separators and on separatrix surfaces is 
likely to be important for solar atmospheric heating. As already mentioned, several 
thousand null points and separators have been found in a single PFSS extrapolation 


made from high-resolution SOLIS data. Using HMI data even more nulls (Edwards 


& Parnell 2015) and associated separators are found . At even higher resolution, we 


are likely to hnd millions of nulls and separators in the solar atmosphere at any one 
time. Additionally, due to the highly-dynamic nature of the solar atmosphere, there 
will be many more topological features than those identified in potential or (non- 
)linear force-free magnetic fields. However, the majority of these additional nulls 
will be located in quiet-Sun regions where the surface magnetic field is immensely 
complex. In active regions, the surface magnetic field is relatively simple, since 
it is dominated by just a few large-scale sunspots or pores. Thus, the magnetic 
field above active regions contains considerably fewer topological features than the 
equivalent sized region of quiet Sun. 
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From our current level of understanding, direct heating at topological features is 
unlikely to be important in active regions, the warmest coronal regions on the Sun. 
Instead, within large (active-region) scale volumes of the solar atmosphere possible 
heating mechanisms (in no particular order) are (i) reconnection associated with the 


complex driving of simple fields (e.g. Galsgaard & Nordlund, 1996 Bowness et al. 
2013 Wilmot-Smith et al. 2010[ Pontin et al.\ ^11|), (ii) reconnection triggered 


by local instabilities in, for example, twisted flux tubes (e.g. Browning et al. [2008 


Bareford et al. 2010 20131 and (iii) the dissipation of waves via phase mixing (e.g. 


Heyvaerts & Priest 19831, resonant absorption (e.g. lonson, 1978) or mode coupling 


(for a review of wave dissipation mechanisms, see De Moortel & Nakariakov 2012). 

In open-field regions, at present it is very difficult to get reliable measurements 
of the magnetic fields. Initial results suggest that isolated coronal nulls that are not 


connected to any separators are fairly common-place in open-field regions (Edwards 


& Parnell, 2015). Until there are more reliable polar-field measurements, the num¬ 


bers of coronal nulls and associated features will remain largely unknown. Thus, it 
is not possible to speculate on whether reconnection at topological features plays a 
major role in the acceleration of the fast solar wind. 

In order to produce a wind that emanates out over the entire area of an open-field 
region, there would have to be an extremely large number of coronal nulls: indeed 
most probably an unfeasably large number. Therefore wave dissipation without null 
points remains the most likely heating/acceleration mechanism in coronal holes (for 
a review, see 


Ofman 2010). 


The topology of the quiet-Sun magnetic field is extremely complex. As the 
resolution and sensitivity of magnetographs increases so too will levels of observable 
mixed polarity held. This will not continue indehnitely and, although not reached 
yet, there is likely to be a limit at which increased resolution does not lead to 
additional small-scale features. So it is possible (depending upon the size, which 
is as yet unknown, of the hnest scale features) that there are sufficient numbers 
of topological features to enable heat due to reconnection and wave dissipation 
mechanisms at topological features to spread over a sufficient area to explain the 
(apparently uniform) heating of the so-called background corona. The question then 
arises: can enough energy be dissipated, by reconnection or waves, at these features? 
Considerably more investigation and theoretical modelling needs to be undertaken 
to answer this. At the moment, numerical models cannot replicate the appropriate 
plasma parameters found on the Sun, so determining the correct energies is difficult. 

The magnetic skeleton associated with the finest-scale photospheric features is 
likely to be on a scale much larger than that found, for instance, in the case where 
(through dynamic reconnection) the local structures about a single separator cas¬ 


cade to small-scales producing many hundreds of separators (e.g. Parnell et al. 


20I0&). In such dynamical cascades that approach/reach a turbulent state, the 


magnetic skeleton may become overly complex and evolve too rapidly for detailed 
knowledge of each and every null point and separator to be of use. However, knowl¬ 
edge of the skeleton before and after this turbulent behaviour could help to identify 
the location of where the turbulent reconnection will occur and the consequences 
of the resulting heating. 

The energy built up over even a relatively short period of slow stressing of the 
field (shearing) and conversion to heat by rapid reconnection (as a result of micro¬ 
instabilities) is likely to be greater than any heating due to wave pulses dissipated 
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at topological features. Furthermore, since waves are trapped only at nulls and 
not at other topological features, the spread of energy through wave dissipation at 
coronal nulls is unlikely to be as great as that from reconnection at coronal nulls, 
separators and separatrix surfaces combined. 

It is worth noting here that, although reconnection leads to localised direct 
Ohmic dissipation, it also launches waves. The energy partitioning of the magnetic 
energy released during reconnection is still a major unanswered question. As already 
mentioned, this question has been tackled from an MHD point of view to ascertain 
whether Ohmic or viscous heating dominates (Longcope & Priest 2007[|Birn et al. 


2009 |Fuentes-Fernandez et al. 2012a[ [Longcope fc Tarr['2012[ Fuentes-Fernandez| 


et al. 20126 Stevenson et al. 20151; there is currently no definitive answer to this 


question. Nonetheless, it is clear that waves can drive reconnection and reconnection 
can launch waves and, therefore, the discussion as to whether waves or reconnection 
is the key heating mechanism is difficult to disentangle, as these mechanisms are 
completely interlinked. 

The question regarding the importance of topological features as the major heat¬ 
ing sites in the solar atmosphere very much depends on the nature of the magnetic 
field within a given region. What the topological features (such as nulls, separa¬ 
tors and, in particular, separatrix surfaces) do is reveal the likely nature of the 
field. Using global models, the boundaries of open-field regions can be identified. 
Furthermore, the volumes of closed field within the open-field regions can be iden¬ 
tified, enabling a much better understanding of the true expanse of the open-field 
volume and also the nature of the magnetic fields within the volume. Both the loca¬ 
tion and geometry of the boundaries of large-scale closed field regions (such as those 
associated with active regions) can also be determined, as can the boundaries of 
the multitude of small-scale closed-field regions that make up the quiet Sun. Thus, 
in terms of coronal heating, knowledge of the magnetic topology is crucial in order 
to identify the nature of the solar atmosphere where different atmospheric heating 
mechanisms are likely to dominate. 

In this article, we have simply considered whether topological features are im¬ 
portant for coronal heating. Thus, we have not considered whether other features, 
such as QSLs, which are principally geometrical in nature, are important for coronal 
heating. It is certainly clear that reconnection can occur in the absence of nulls and 


separators (e.g. 

Schindler et al. 

1988 

Hesse & Schindler 1988 

Galsgaard & Nord- 

lund 

[19961 [Aulanier et al.\ 2005 2006 Wilmot-Smith et al.\ 

2010 Pontin et al. 

2011 

1, as can the dissipation of waves (e.g. lonson 1978 Heyvaerts & Priest 1983). 


Furthermore, we have not addressed the question of whether topological features 
are important for other solar events, such as solar flares, CMEs, X-ray jets or X-ray 
bright points. There are at present a number of published works which strongly 


suggest that they are indeed important (e.g. 

Parnell et al. 

1994 

Aulanier et al. 

2000 

Ugarte-Urra et al. 2007 

Pariat et al. 2009 

Masson et al. 

20091. 
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